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Motivated by recent experiments on the quantum-spin-liquid candidate material LiZn2Mo308, 
we study a single-band extended Hubbard model on an anisotropic Kagome lattice with the 1/6 
electron hlling. Due to the partial hlling of the lattice, the inter-site repulsive interaction is necessary 
to generate Mott insulators, where electrons are localized in clusters, rather than at lattice sites. 

It is shown that these cluster Mott insulators are generally U(l) quantum spin liquids with spinon 
Fermi surfaces. The nature of charge excitations in cluster Mott insulators can be quite different 
from conventional Mott insulator and we show that there exists a novel cluster Mott insulator where 
charge fluctuations around the hexagonal cluster induce a plaquette charge order (PCO). The spinon 
excitation spectrum in this spin-liquid cluster Mott insulator is reconstructed due to the PCO so 
that only 1/3 of the total spinon excitations are magnetically active. Based on these results, we 
propose that the two Curie-Weiss regimes of the spin susceptibility in LiZn 2 Mo 3 08 may be explained 
by hnite-temperature properties of the cluster Mott insulator with the PCO as well as fractionalized 
spinon excitations. Existing and possible future experiments on LiZn 2 Mo 3 08, and other Mo-based 
cluster magnets are discussed in light of these theoretical predictions. 


I. INTRODUCTION 

If there is no spontaneous symmetry breaking, the 
ground state of a Mott insulator with odd number of 
electrons per unit cell may be a quantum spin liquid 
(QSL)P1. QSL is an exotic quantum phase of matter with 
a long-range quantum entanglement^ and is character¬ 
ized by fractionalized spin excitations and an emergent 
gauge structures at low energie^^^^. It has been suggested 
that some frustrated Mott insulating systems which are 
proximate to Mott transitions may provide physical re¬ 
alizations of These QSLs are believed to arise 

from strong charge fluctuations in the weak Mott regime, 
which can generate sizable long-range spin exchanges or 
spin ring exchanges and suppress possible magnetic or- 
derings^. Several QSL candidate materials, such as the 
2D triangular lattice organic materials /t:-(ET) 2 Cu 2 (CN )3 
and EtMe 3 S b[Pd(d mit) 2 ] 2 r and a 3D hyperkagome sys¬ 
tem Na 4 lr 30 jmti^, are expected to be in this weak Mott 
regime. These weak Mott-insulator U(l) QSLs are ob¬ 
tained as a deconfined phase of an emergent U(l) lattice 
gauge theor}®^, where the electron is fractionalized into 
spin-carrying spinons and charged bosons. The charge 
excitations are gapped and the low-energy physics of the 
QSLs is described by a spinon Fermi surface coupled to 
the emergent U(l) gauge field. 

Recently LiZn2Mo308 has been proposed as a new 
candidate material for QSL ground statd^^^hJO^ Besides 
the usual QSL phenomenolog}ff^®2l^ the experiments re¬ 
veal two Curie-Weiss regimes of the spin susceptibility at 
finite temperatures. The low temperature Curie-Weiss 
regime is governed by a much smaller Curie-Weiss tem¬ 
perature than the high temperature one and a reduced 



FIG. 1. (Color online.) (a) The electron configuration in the 
CMI without the PCO when V 2 ^ ti and Vi ^ t 2 - (b) The 
electron conhguration in the CMI with the PCO. Three elec¬ 
trons hop resonantly in each hexagon that is marked by a 
(red) circle. These marked hexagons form an emergent trian¬ 
gular lattice (with lattice vectors ai,a 2 ). 


Curie constant which is 1/3 of the high temperature 
counterpart. In this work, to understand how one can 
achieve the QSL phenomenology and the puzzling two 
Curie-Weiss regimes in this material, we consider a 1/6- 
fllled extended Hubbard model with nearest-neighbor re¬ 
pulsions on an anisotropic Kagome lattice. We first pro¬ 
pose the existence of plaquette charge order (PCO) for 
the charge degree of freedom. We emphasize that the 
emergence of the PCO reconstructs the spin state of the 
system. We further propose a U(l) QSL with spinon 
Fermi surfaces for the spin ground state and a PCO for 
the charge ground state in this system. The Mott in¬ 
sulators in partially filled s ystem s arise due to the large 
nearest-neighbor repulsionsP^ and localization of the 
charge degrees of freedom in certain cluster units. Hence 
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such Mott insulators are dubbed “cluster Mott insula- 
tors” (CMIs;PHlIl. 

The single-band extended Hubbard model is defined 
on the anisotropic Kagorne lattice of the Mo sites (see 
Fig. and is given b}® 

H= ^ Hi(ctc^v + ^-c-) +Vnin^] 

{ij)eu 

+ ^ H2(4<TS<T + ^-'^-) + ^2ninj] 

{ij)ed 

+ '2^^^ ~ 2 ^^’ 

i 

where the spin index a is implicitly summed, ct^ 

(c-^) creates (annihilates) an electron with spin a at lat¬ 
tice site and ti, Vi and ^ 2 , V 2 are the nearest-neighbor 
electron hopping and interaction in the up-pointing tri¬ 
angles (denoted as ‘u’) and the down-pointing triangles 
(denoted as ‘d’), respectively, rii = elec¬ 

tron occupation number at site i. Since there exists only 
one unpaired electron in each Kagorne lattice unit celF^, 
the electron filling for this Hubbard model is 1/6. 

Although the down-triangles are larger in size than the 
up-triangles in LiZn 2 Mo 308 , because of the large spatial 
extension of the 4d Mo electron orbitals we think it is 
necessary to include the inter-site repulsion V 2 for the 
down-triangles. For LiZn2Mo308 we expect ti > t 2 and 
U > Vi ^ K 2 - Because of the very dilute electron fill¬ 
ing, although the Hubbard U is the largest energy scale, 
it alone can only remove double electron occupation on 
a single lattice site and cannot localize the electron. If 
there is no Vi or V 2 interactions, the electrons can still 
transport on the lattice without encountering any elec¬ 
tron double occupancy on a single lattice site. So we 
need Vi and V 2 to localize the electrons in the (elemen¬ 
tary) triangles of the Kagorne lattice instead of the lattice 
sites. 

Let us first explain the electron localization in the 
absence of V 2 - Clearly, as t 2 is the hopping between 
the up-triangles, when Vi ^ ^ 2 , the electrons are local¬ 
ized on the up-triangles with one electron per up-triangle 
(see Fig. Ef)- In this picture, the localized electron can 
hop freely among the three lattice sites within each up- 
triangle and gain local kinetic energy ^ 0{ti) while the 
electron number on the down-triangle is strongly fluctu¬ 
ating. After V 2 is introduced, as V 2 increases, the con¬ 
figuration with more than one electrons on the down- 
triangles (like the one in Fig. [^) becomes less favor¬ 
able energetically. When the interaction energy cost 
(^ 0 (V 2 )) on the down-triangle overcomes the local ki¬ 
netic energy gain ^ 0(ti), the electron number on each 
down-triangle is also fixed to one and the electrons can 
no longer move freely within each up-triangle. Instead, 
the electrons develop a collective motion. For instance, 
in Fig.[^, the three electrons on the hexagon at position 
R can tunnel between the configuration occupying sites 
1,3,5 and the other configuration occupying sites 2,4,6. 


This collective electron tunnelling process preserves the 
electron number on each triangle and is the dominant 
physical process below the Mott gap. We show this col¬ 
lective electron tunnelling gives rise to a long-range PCO 
that breaks the lattice symmetry spontaneously. With 
the PCO, the electrons are preferentially tunnelling back 
and forth on the hexagons that are marked with a (red) 
circle (see Fig. Ef)- We will refer these special hexagons 
as “resonating” hexagons. On these resonating hexagons, 
the three electrons form a linear superposition state of the 
two electron configurations with sites 1,3,5 or sites 2,4,6 
occupied. We emphasize and will show in Sec. [H] that 
the emergence of the PCO in the CMI is a quantum ef¬ 
fect and cannot be obtained from the classical treatment 
of the electron interaction. 


With the PCO, 1/3 of the elementary hexagons be¬ 
come resonating. As shown in Fig. [^, these resonating 
hexagons form an emergent triangular lattice (ETL). The 
PCO triples the original unit cell of the Kagorne lattice, 
and the localized electron number in the enlarged unit 
cell now becomes 3, which is odd. According to Hastings’ 
theorenJD^ the CMI with the PCO is not connected to a 
trivial band insulator and the QSL is expected. In the 
resulting U(l) QSL, we obtain 9 mean-field spinon sub¬ 
bands, compared to the 3 spinon bands in the U(l) QSL 
for the CMI without the PCO. The 9 mean-field spinon 
sub-bands are obtained by splitting the 3 spinon bands of 
the CMI without the PCO, and this is the reason why we 
use the term ^^sub-bands’^ A direct band gap separates 
the lowest spinon sub-band from other spinon sub-band 
in the presence of the PCO. The lowest spinon sub-band 
is completely filled by 2/3 of the spinons, leaving the re¬ 
maining 1/3 of the spinons to partially fill the second and 
third lowest spinon sub-band. Because of the band gap, 
the only active degrees of freedom at low energies are 
the spinons in the partially filled spinon sub-band, and 
the fully-filled lowest spinon sub-band is inert to external 
magnetic field at low temperatures as long as the PCO 
persists. Therefore, only 1/3 of the magnetic degrees of 
freedom are active at low temperatures. If one then con¬ 
siders the local moment formation starting from the band 
filling picture of the spinons (just like electrons occupy¬ 
ing the same band structure) only the 1/3 of the spinons 
from the partially filled upper bands would participate 
in the local moment formation. This means the CMI 
with the PCO would be continuously connected to the 
Curie-Weiss regime with the 1/3 Curie constant at finite 
temperature (compared to the case when all spinons can 
participate in the local moment formation). This would 
explain the two Curie-Weiss regimes in the spin suscep¬ 
tibility data of LiZn2Mo308. More precise connection to 
the existing experiments is discussed later. 


The rest of the paper is structured as follows. In 
Sec.|llj we show the CMI develops the PCO in the charge 
sector when every triangle contains only one electron. 
We generalize the Levin-Wen variational string wave- 
function approaclP^ to study the reconstruction of the 


spinon band structure by the PCO in Sec. HI In Sec. IV 
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we explain the consequence of this reconstructed spinon 
band structure and discuss the low-temperature mag¬ 
netic susceptibility. In Sec. |V| we connect our theory to 
the experiments on LiZn 2 Mo 308 , suggest possible future 
experiments, and discuss other Mo based cluster mag¬ 
nets. Finally, we discuss the quantum chemistry justifi¬ 
cation of the extended Hubbard model in Appendix A 
complementary explanation of the double Curie regimes 
based on the spin state reconstruction for the strong PCO 
regime is given in Appendix In Appendix we pro¬ 
vide the detailed formalism of the mean-field theory in 
the main text. 


II. THE EMERGENCE OF THE PLAQUETTE 
CHARGE ORDER 

As a preparation step, we first employ the standard 
slave-rotor representation and map the low-energy charge 


sector Hamiltonian into a quantum dimer model on the 
dual honeycomb lattice. With the quantum dimer model, 
we then show the system should develop the PCO. 


Using standard slave-rotor represent at iorP^, we first 
express the electron operator, cj^ = where the 

bosonic rotor carries the electron charge and the 

fermionic spinon (//^) carries the spin quantum num¬ 
ber. As it is well-known, the slave-rotor representation 
enlarges the Hilbert space. To constrain the enlarged 
Hilbert space, we introduce an angular momentum vari- 
able L|, L| = [E<t A/ial “1/2, where L| is conjugate 
to the rotor variable with = i6ij. In terms of 

the slave-rotor variables, the extended Hubbard model is 
now reformulated as 


{ij)eu 

+ E + h.c.) + V 2 {Lt + b(q + b] 

{i3)ed 

+ E + 2 ~ E A/icr)]; (2) 

i (J 

where we have introduced hi as a Lagrange multiplier to impose the Hilbert space constraint, since the on-site 
interaction U is assumed to be the largest energy scale, in the large U limit the double electron occupation is always 
suppressed. Hence, the angular momentum variable Lf primarily takes Lf = 1/2 ( — 1/2) for a singly-occupied (empty) 
site. 

Via a decoupling of the electron hopping terms in H into the spinon and rotor sectors, we obtain the following two 
coupled Hamiltonians for the spin and charge sectors, respectively, 

Hs = - + h.c.) 

{ij) * 

Fc = E [-2J^cos(0i - 0^) + + i)(q + ^)] + E + ^)] ’ 

(*i> * 

where J?f = and Uj = h (^2), Vij = V (V2) for the 

bond ij on the up-triangles (down-triangles). Here, we have chosen the mean-field ansatz to respect the symmetries 
of the Kagome lattice. The Hamiltonians Hg and He are invariant under an internal U(l) gauge transformation, 
fia ^ud ^ aij + Xi “ Xj- 

We now focus on the charge sector and study the low energy physics of the charge sector. From the previous 
slave-rotor formulation, the charge sector Hamiltonian is given by 

Hc = Y. -2-/^cos(d, - 0^) + + i)(q + ^) + E him + (5) 

hj) « 

where we have dropped the U interaction term with the understanding that Li = ±1/2 in the large U limit. This 
charge sector Hamiltonian can be thought as a Kagome lattice spin-1/2 XXZ model in the presence of an external 
magnetic field upon identifying the rotor operators as the effective spin ladder operators, = Lf where Lf\Ll = 
=f|) = \Li = ±^). Thus the corresponding effective spin-L model reads 

= E + h.c.) +if q] + ft ^ if, 

{ij) i 


( 3 ) 

( 4 ) 


( 6 ) 
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FIG. 2. (Color online.) The two collective hopping processes that contribute to the ring electron hopping or the ring exchange 
in Eq. Q. The (blue) solid ball represents the electron or the charge rotor. 


in which we have made a uniform mean-field approximation such that + 3(Vi + V 2 ) = h. The 1/6 electron filling is 
equivalent to Ns~^ M where Ng is the total number of Kagome lattice sites. 



FIG. 3. (Color online.) The anisotropic Kagome lattice and 
the dual honeycomb lattice (in gray). 

As we explained in Sec.jTj when Vi {V 2 ) is large enough 
compared to ^2 (ti), the electron number on each trian¬ 
gle, both up-triangle and down-triangle, is fixed to be 
one. The electron occupation configuration that satisfies 
this condition is highly degenerate. The presence of the 
electron hopping, i.e. LJ , lifts this classical degener¬ 
acy and the effective interaction can be obtained from a 
third-order degenerate perturbation theory. The result¬ 
ing effective ring exchange Hamiltonian is given as 

^c,ring = ~ E Jring(E^Z /2 Lq + h.c.)^ (7) 

O 

where “o” refers to the elementary hexagon of the 
Kagome lattice, Jring = /V 2 + /V^ and 

“1, • • •, 6” are the 6 vertices on the perimeter of the ele¬ 
mentary hexagon (see Fig. |^. This ring Hamiltonian in 
Eq. describes the collective tunnelling of three electron 
charges between A and B configurations in Fig. 

We now map i7c,ring into a quantum dimer model on 
the dual honeycomb lattice that is formed by the centers 
of the triangles on the Kagome lattice (see Fig. [^. As 


depicted in Fig.[^ a dimer is placed on the corresponding 
link if the center of the link (or the Kagome lattice site) 
is occupied by an electron charge. The rotor operator 
Lf simply adds or removes the charge dimer. So i7c,ring 
is mapped into the quantum dimer model with only a 
resonant term, 

-f^c,ring — J ring ( | Oa ) (Ob | T |Ob)(Oa|) ( 8 ) 

0 

where |Oa) and |Ob) refer to the two charge dimer cover¬ 
ing configurations in the elementary hexagon “O” of the 
dual honeycomb lattice as shown in Fig. 

In Ref.[23l Moessner, Sondhi and Chandra studied the 
phase diagram of the quantum dimer model on the hon¬ 
eycomb lattice quite extensively. In the case with only 
the resonant term of our model in Eq. (|^, they found 
a translational symmetry breaking phase with a plaque- 
tte dimer order, in which the system preferentially gains 
dimer resonating (or kinetic) energy through the resonat¬ 
ing hexagons on the dual honeycomb lattice (see Fig.|^). 
The dimers on resonating hexagons form a linear super¬ 
position of the dimer covering configurations |Oa) and 
I Ob). This indicates that our model is unstable to trans¬ 
lational symmetry breaking. The plaquette dimer order 
of the quantum dimer model is then mapped back to 
the plaquette charge order (PCO) on the Kagome lat¬ 
tice (see Fig. &)■ Just like the resonating dimers, the 
three electron charges on the resonating hexagons also 
form a linear superposition of two occupation configura¬ 
tions in Fig|^. This is a quantum mechanical effect and 
cannot be obtained by treating the inter-site electron in¬ 
teractions Vi and V 2 in a classical fashion. Moreover, this 
PCO can be regarded as a local charge resonating valence 
bond (RVB) stat e whi ch contrasts with the spin singlet 
RVB of Anderso] j^^F5 [ 
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(a) 



works. 

With the PCO, the electrons are preferentially hopping 
around the perimeters of the resonating hexagons on the 
Kagome lattice. These resonating hexagons are period¬ 
ically arranged, forming an emergent triangular lattice 
(see Fig. Ef)- Due to the translational symmetry break¬ 
ing, this emergent triangular lattice (ETL) has an en¬ 
larged unit cell that includes 9 sublattices compared to 
3 sublattices in a Kagome lattice (see Fig. [^. 


FIG. 4. (Color online) (a) The plaquette charge dimer or¬ 
dering pattern on the dual honeycomb lattice. The charge 
dimers have a high probability to occupy the bold bonds of 
the resonating hexagons (see the main context), (b) The cor¬ 
responding PCO on the Kagome lattice. We mark the res¬ 
onating hexagons with both dark bonds and the red circles. 

We note that similar type of PCO has already been 
obtained for extended Hubbard models with fermions or 
hard-core bosons on an isotropic Kagome lattice with 1/3 
and 2/3 fillings in certain parameter regimes in previous 
works^^H^. The result was obtained either through per- 
turbatively mapping to the quantum dimer model or by a 
Hatree-Fock mean-field calculation. In particular, Ref.fSOl 
applied the quantum Monte Carlo technique to simulate a 
hardcore boson Hubbard model on an isotropic Kagome 
lattice and discovered the PCO for 1/3 and 2/3 boson 
fillings. Because our model is defined on the anisotropic 
kagome lattice, it is not exactly the same as the previous 


III. THE SPINON BAND STRUCTURE 

In Sec. [Hj using a slave-rotor approximation, we have 
shown that the system is unstable to the development of 
the PCO in the CMI where both up-triangle and down- 
triangle contain only one electron. Although our result is 
obtained by first starting from a translationally invariant 
mean-field ansatz, as we argue below, the PCO breaks 
the translational symmetry of the spinon mean-field state 
and the modified spinon band structure makes the PCO 
even more stable^. 

A. Spin charge coupling in the CMI with the PCO 

To understand how the PCO in the charge sector influ¬ 
ences the spinon sector, we first consider the low-energy 
effective ring hopping model in the CMI where both up- 
triangle and down-triangle contain only one electron. 


Hring = - E E ['*^l(4aC6a4/3C4/}47% + . (9) 

O a/37 


where Ki = 6 tf/k 2 ^ and IK 2 = are readily obtained from the third-order degenerate perturbation theory. Here, 

are the 6 vertices in the elementary hexagon of the Kagome lattice. 

Using the slave-rotor representation in Sec. for the electron operator cj^ = the ring hopping 

model i^ring can be decoupled as 

Hring = L 2 L'^Lq X Miq^ 4S2 h.C.) -\-K2{L^L 2 L'^Lq X Mi23456 + ^-C.)] (10) 

o 

= [lKi( I Oa)(Ob 1^165432 + I Ob )(Oa 1^165432) + ^^2 (|Oa) (Ob 1^123456 + |Ob) (Oa 1^123456)] : (H) 

0 

where |Oa) and |Ob) are the two charge dimer coverings in Fig. Here we are focusing on the charge sector and 
treating the spinon sector in a mean-field fashion, i.e. 


^ijklmn — E ifLfjJkpfllsfll'lfrn), 


a/37 


( 12 ) 


where the lattice sites are arranged either clockwise or anti-clockwise and is the complex con¬ 

jugate of Mijkimn- What we did in Eq. 0 is to directly couple the charge sector quantum dimer model with the 
spinon sector. By doing this, we can study how the charge sector is influenced by the spinon sector. 

The parameter Mij^imn is evaluated in the spinon mean-field ground state, which we explain below. Eor a time- 
reversal invariant system, we expect Mijkimn = ^ijkimn' assume the spinon sector respects the translational 

invariance of the Kagome lattice, the resulting charge sector model would be equivalent to i^c,ring iii Eq. 0 and also 
to the quantum dimer model in Eq. 0 except for the renormalized couplings, and this would immediately imply 
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FIG. 5. (Color online) (a) The Kagome lattice is partitioned into unit cells on the ETL. The unit cell (marked by a gray 
triangle) contains 9 sublattices that are labelled by “1,2,- • • ,8,9”. (b) The spinon hoppings on the bonds that surround around 
a resonating hexagon. 


the system should develop the PCO and spontaneously break the translational symmetry of the Kagome lattice. In 
turn, the breaking of lattice symmetry by the PCO would then influence the spinon sector. To understand how the 
spinon band structure is modified by the underlying PCO in the charge sector and how the modified spinon band 
structure feeds back to the charge sector, we take the enlarged unit cell of the ETL and introduce the following spinon 
mean-field (hopping) Hamiltonian (see Fig. If)- 


— E] [ *i(/r1ct/r6ct + /L<t/r3<t + /L(t/r5(t) *2 (/L<t/r2(t + /L<t/r4<t + /L<t/r6<t) 


R 


(/rIct/rTct 

«^R9cr/R-|-ai ,2cr «^R7cr«^ R-|-a2,3cr 

where R labels the unit cell of the ETL and “1, 2 


+ /rGct/rTct + // 
ft 


R2cr/R8cr + /rSct/rSct + fil9afll4a + fil9afll5a) ^2 (// 


t/fI 


/R7cr«^R-|-a2,4cr 


+ /rSct/r —ai—a2,5cr 


+ /rSct/r — ai —a 2 , 6 cr 


RQcr-' R-|-ai,lcr 

h.c.]. 


(13) 


•, 8, 9” label the 9 sublattices of the ETL. This choice of spinon 
hopping parameters respects the 3-fold rotation symmetry and the reflection symmetry of the resonating hexagons 
(see Fig. [^). Moreover, in Eq. (13), the spinon hoppings are related to the charge sector via 


=ii(L+(R)L6(R)) 

(14) 

=i2(L+(R)L2-(R)) 

(15) 

i\=ti{Lt{K)LT{K)) 

(16) 

t '2 — '^2{Lq (R-)ii (R + ^i)) 

(17) 


such that the influence of the charge sector on the spinon sector is captured. 


In the presence of the PCO, we expect ti > i'^ and 
^2 > i '2 due to the presence of the PCO. According 
to Eqs. (p^, ([T^, ([ 1 ^ and ([T7|), the presence of the 


PCO would enhance the bonding of the charge rotors and 
then the spinon hoppings in the resonating hexagons and 
weakens the ones in the non-resonating hexagons. The 
enhanced spinon hoppings in the resonating hexagons 
further strengthen the couplings of the ^ring in the res¬ 
onating hexagons through Mij^imn in Eq. (11). Thus the 


PCO would become more stable if the coupling between 
spinon and charge excitations is switched on. 


B. Generalized Levin-Wen’s variational dimer 
wavefunction approach and the spinon band 
structure in the presence of the PCO 

We now consider the combination of the spinon hop¬ 
ping model in Eq. (13) with the ring hopping model 


iLring in Eq. ( |lT] ). It was pointed by Levin and WerP^^^ 
that quantum dimer model is an example of string-net 
models. In Ref. [22l Levin and Wen developed a vari¬ 
ational string wavefunction approach (or string mean- 
field theory) to describe the properties of quantum dimer 
model. To make the nomenclature consistent, we refer 
Levin-Wen’s variational string wavefunction approach as 
variational dimer wavefunction approach in the follow¬ 
ing. Since the charge sector is described by a quan¬ 
tum dimer model, we can extend Levin-Wen’s varia¬ 
tional string wavefunction approaclP^ to solve the cou- 
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pled charge and spinon problem in Sec III A In Levin and 
Wen’s original work, the variational dimer wavefunction 
approach was designed for pure quantum dimer model. 
The new ingredients of our problem are the presence of 
the spinon degrees of freedom, and the coupling and the 
mutual feedback between the spinous and charge dimers. 

We describe below the variational dimer wavefunction 
approach that is used to optimize the Hamiltonian Hring 
in Eq. (11) for the charge dimers. Following Levin and 
Wen, the variational dimer wavefunction is parametrized 
by a set of variational parameters {zi} where Zi is defined 
on each link of the dual honeycomb lattice. Here the links 
on the dual honeycomb lattice are also parametrized by 
the Kagome lattice sites that are located at the centers of 
the links. These variational parameters Zi are also termed 
as string (or dimer) fugacity by Levin and WerP^. For 
each set of { 2 :^}, the variational dimer wavefunction is 
given by 




TT |Q)i + Zi\l)i 


(18) 


where |0)i and |l)i define the absence and presence of the 
electron charge at the Kagome lattice site i or the dimer 
on the corresponding link on the dual honeycomb lattice, 
respectively. Moreover, we have the following relations 
by definition. 


the Hamiltonian iLc,ring with respect to T({zi}) and op¬ 
timizing the energy by varying z. This static variational 
approach, however, cannot directly produce the plaque- 
tte ordered phase of the quantum dimer model. What 
it gives is a translationally invariant variational ground 
state. To obtain the right result. Levin and Wen de¬ 
veloped a dynamical variational approach. Namely, for 
each static variational ground state, one checks the sta¬ 
bility of the variational phase by considering the quan¬ 
tum fluctuation of this phase. In the model that Levin 
and Wen were considering^, they found some modes in 
a translationally invariant variational ground state can 
become unstable and drive a dimer crystal ordering. We 
expect similar physics should happen to our quantum 
dimer model iLc,ring- 

Unfortunately, the dynamical variational approach by 
Levin and Wen is not a self-consistent variational ap¬ 
proach and cannot be extended to the combined spinon 
and charge dimer problem that we are interested in here. 
Since we know our quantum dimer model i7c,ring gives 
the ground state with the PCO and we have argued that 
coupling the charge (dimer) with the spinous makes the 
PCO even more stable in Sec. we now introduce 

the PCO into the system by explicitly breaking the lat¬ 
tice symmetry. That is, we modify the ring hoppings Ki 
and IK 2 in ^ring of Eq. (10). For the resonating hexagons, 
we change 


ni|0)i = 0, ni\l)i = \l)i, (19) 

L+| 0 ), = | 1 ),, Lr|l), = | 0 )„ ( 20 ) 

where rii counts the electron number (or the number of 
dimers) at site i. 

We employ the symmetry of the PCO to reduce the 
number of free variational parameters in the dimer wave- 
function ^({ 2 ;^}). Using the symmetries of the ETL, we 
find that only two variational parameters are needed 

^i(R) = 2^2 (R) = 2:3 (R) = Z4{'R) 

= z^CR) = 2;6(R) = 2 ;, (21) 

2;7(R) = 2;8(R) = 2;9(R) = z, (22) 

where 2 ;^(R) (/i = 1,2,-•• ,9) refers to the variational 
parameter of the fith sublattice at the unit cell R (see 
Fig. §. We have reduced the set of variational param¬ 
eters in the variational dimer wavefunction to 2 ; and 5. 
Moreover, 2 ; and z are not independent from each other. 
This is because of the charge localization constraint, i.e. 
every triangle contains only one electron. In terms of the 
dimer language, this constraint is that every dual hon¬ 
eycomb lattice site is connected by only one dimer. To 
satisfy this constraint, we only need to require 

(ni(R)) + (n6(R)) + (n 7 (R)) = 1, (23) 

where the expectation value is taken for the variational 
wavefunction ^({ 2 :^}). This relation connects z to z. 

For the quantum dimer model i7c,ring in Eq. (|^, vari¬ 
ational (or mean-field) phase is obtained by evaluating 


Ki IKi(l + IK2 1^2(1 + ( 5 ), ( 24 ) 

and for the non-resonating hexagons, we use 

IKi — y IKi(l — (i), IK2 — y 1^2(1 — (^), 

where 6 (with > 0) is a phenomenological parame¬ 
ter that breaks an appropriate lattice symmetry for the 
PCO. This modification of the ring hoppings captures 
the spatial modulation of the energy in the system when 
the PCO is present. This phenomenological way of in- 
trodncing the PCO is very similar in spirit to Henley’s 
approaclP^to the order by disorder, where a phenomeno¬ 
logical interaction is introduced into the energy or the 
free energy to model the ground state selection due to 
the quantum fluctuation. 

We now solve the combined Hamiltonian of and 
^ring with the modified ring hoppings self-consistently. 
This self-consistent approach is expected to underesti¬ 
mate the PCO for any fixed S and thus underestimates 
the reconstruction of the spinon band structure due to 
the PCO. Nevertheless, to understand the generic fea¬ 
tures of the spinon band structure reconstruction in the 
presence of the PCO, we can simply vary the phenomeno¬ 
logical parameter 6 and study the spinon band structure 
from this self-consistent approach. 

The evolution of the mean-field spinon band structure 
is depicted in Fig. When ^ = 0, there is no PCO 
and the symmetry of the Kagome lattice is preserved. 
The spinon band structure contains 3 bands in the 1st 
Brioullin zone of the Kagome lattice (BZl in Fig. [^). 
The 3 spinon bands are well separated in energy and 









(a) (5 = 0 (b) S = 0.3 




(d) Brioullin zones 


FIG. 6. (Color online.) (a-c) The evolution of the spinon 
sub-bands as S is varied. In the figure, we choose Ki = 
4, IK 2 = 1 and ti = l,t 2 = 0.5. The (red) dashed line is 
the Fermi energy of the spinons. The energy unit is set to 
ti. The (red) arrow indicates a 2-fold degeneracy. The de¬ 
tails of the band structures are discussed in the text, (d) 
The large and small hexagons define the Brioullin zones 
of the Kagome lattice (BZl) and the ETL (BZ2), respec¬ 
tively. Setting the Kagome lattice constant to unity, we have 
K = (f = (f and qi = (f ,0),q. = 


3 
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have no direct nor indirect overlap, and we can simply 
focus on the lowest band as the spinons only fill half 
of the lowest band. So in Fig. we only need to plot 
the evolution of the lowest spinon band. In Fig. [^, we 
have further folded the lowest spinon band onto the 1st 
Brioullin zone of the ETL (BZ2 in Fig. [^) and obtain 
3 spinon sub-bands. We use '‘’‘spinon sub-bands’^ to refer 
the spinon bands plotted in the BZ2 of the ETL. 

In the presence of the PCO for a finite 6, the system 
has 9 sublattices (see Fig.[^, the 3 spinon bands at ^ = 0 
are further split into 9 spinon sub-bands, and the lowest 
spinon band at (5 = 0 is split into 3 spinon sub-bands. 

We now explain the evolution of the spinon bands as 
the PCO is enhanced by increasing the variational pa¬ 
rameter 6. At 6 = 0, there is no PCO, and the 2nd and 
the 3rd spinon sub-bands touch at the zone boundary 
of the BZ2. A finite 6 creates the PCO and breaks the 
translational symmetry of the Kagome lattice. We have 
ii > i'l and ^2 > as previously expected. The band 
touching of the 2nd and the 3rd spinon sub-bands at the 
zone boundary of the BZ2 is lifted by the level repulsion 
(see Fig. IT)- A direct band gap is created between the 
lowest spinon sub-band and upper spinon sub-bands. We 
emphasize this feature is generic and is not specific to the 
ring hoppings and electron hoppings that are chosen in 
Fig. and we have also explicitly checked many other pa¬ 


rameter choices. As the parameter 6 is further increased 
and the PCO becomes even stronger, the direct band gap 
gets larger and larger, and eventually the lowest spinon 
sub-band is fully separated from the other sub-bands by 
a full band gap (see Fig. |^). Therefore, the band gap of 
the spinons is set by the stiffness of the PCO. 

Even though the PCO enlarges the unit cell from 3 sites 
of the Kagome lattice to 9 sites of the ETL, the spinon 
Fermi surface always exists. This is because the number 
of electrons or spinons per unit cell is 3 with the PCO. 
Because of the direct band gap, the lowest spinon sub¬ 
band is completely filled by the spinons which comprise 
2/3 of the total spinon number, and the remaining 1/3 of 
spinons partially fill the upper 2 sub-bands and give rise 
to the spinon Fermi surfaces. Therefore, the interna l U(l) 
gauge field is expected to be in the deconfining phas^^^^, 
and we obtain U(l) QSL for the ground state. 


IV. THERMAL TRANSITION AND SPIN 
SUSCEPTIBILITY 

Because the PCO breaks the lattice symmetry, this 
implies that there exists a thermal phase transition at a 
finite temperature which destroys the PCO and restores 
the lattice symmetry. This thermal transition is expected 
to occur at T* ^ O(IKi) = OitXjV^) (because Ki ^ 
IK 2 ) when the local electron resonance in the elementary 
hexagons loses the quantum phase coherence. 

Based on the understanding of the spinon band struc¬ 
ture in Sec. m we describe the behavior of the spin sus¬ 
ceptibility in the low temperature regime with the PCO 
(T < T*) and in the high temperature regime without 
the PCO (T > T*). Since the U(l) QSL ground state 
has spinon Fermi surfaces, we expect a finite (Pauli-like) 
spin susceptibility in the zero temperature limit. At finite 
temperatures, one should recover the Curie-Weiss law for 
the spin susceptibility. What are the Curie constant and 
the Curie-Weiss temperature that characterize the Curie- 
Weiss law for T < T* and T > T*? The Curie constant 
measures the number of the active local moments. Let 
us now consider the local moment formation regime for 
T < T*. As long as the PCO is not destroyed by ther¬ 
mal fluctuations which is the case for T < T*, the direct 
band gap between the lowest spinon sub-band and up¬ 
per sub-bands would be present, and the lowest spinon 
sub-band is fully filled by 2/3 of the spinon numbers. A 
fully-filled spinon is inert to an external magnetic field, 
and thus, only the 1/3 of the spinons from the partially 
filled upper sub-bands contribute to the local moment, 
which comprise 1/3 of the total number of electrons in 
the system. Therefore, the low temperature Curie con¬ 
stant in a DC susceptibility measurement for T < T* 
is 


^ + 1 ) Na 

3kB 3 


(26) 


where g ^ 2 is the Lande factoifl^^, s = 1/2, and Na 
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is the total number of up-triangles in the system. From 
the electron filling fraction, we know equals the to¬ 
tal electron number N^. Because only 1/3 of the total 
spins are responsible for the low-temperature magnetic 
properties, the Curie constant is only 1/3 of the one at 
very high temperatures where all the electron spins are 
supposed to be active. 

For T > T*, the PCO is destroyed by thermal fluctua¬ 
tion, and the direct band gap between the lowest spinon 
sub-band and the upper spinon sub-bands is closed. All 
the localized electrons are active and contribute to the 
local moment, and thus the Curie constant in this high 
temperature regime is 




3/cb 


Na. 


(27) 


which is 3 times the low temperature one, C^. 

As for the Curie-Weiss temperature, it is hard to make 
a quantitative prediction from the spinon Fermi sur¬ 
face. But it is noted that the Curie-Weiss temperature is 
roughly set by the bandwidth of the active spinon bands 
in the QSL phase. At low temperature PCO phase, the 
active spinon bands are the partially-filled upper spinon 
sub-bands on the ETL (see Fig.[^). As one can see from 
Fig. [^, the bandwidth of the active spinon bands is sig¬ 
nificantly reduced when the PCO is present compared to 
Fig. when the PCO is absent. As a result, we expect 
a much reduced Curie-Weiss temperature in the presence 
of the PCO at T < T* compared to the case in the ab¬ 
sence of the PCO at T > T*. In the absence of the PCO 
at T > T*, as all the spinon sub-bands are active, the 
Curie-Weiss temperature is set by the total spinon band¬ 
width in Fig. [^. Finally, we want to point out that the 
double Curie regimes in the spin susceptiblity is a finite 
temperature property and independent from whether the 
spin ground state is a spinon Fermi surface U(l) QSL or 
not. It is the PCO that reconstructs the spin states and 
leads to the double Curie regimes. In the Appendix \B\ 
we provide a complementary explanation of the double 
Curie regimes from the spin state reconstruction in the 
strong PCO regime. 


V. DISCUSSION 
A. Applications to LiZn2Mo308 

As we discuss in Sec. |IV[ the CMI with the PCO 
provides two Curie-Weiss regimes in spin susceptibility. 
Armed with these results, we here propose that the Mo 
system in LiZn2Mo308 may be in the CMI with the PCO 
at low temperatures. Because the PCO triples the unit 
cell, the thermal transition at T* is found to be first or¬ 
der in a Landau theory analysis for a clean systenP^. 
In reality, LiZn2Mo308 is influenced by various disor¬ 
ders or impurities (e.g. the mixed Li/Zn sites and mo¬ 
bile Li ions]P^. For example, impurities would broaden 


the charge ordering transitioiP^. This may explain why 
a sharp transition is not observed in the experiments^. 
Nevertheless, the experiments do observe a peak around 
lOOK in heat capacit^ttSl which might be related to the 
smeared-out phase transition. 

Based on the fact that there is no obvious ordering 
down to ^ O.IK for LiZn2Mo308 and the app arent gap¬ 
less spin excitation in neutron scatteringi^Hi^, we further 
propose that the system is in the U(I) QSL with spinon 
Fermi surfaces (as well as the PCO) of the CMI that is 
obtained in Sec. ^ With the spinon Fermi surfaces, we 
expect the usual behaviors of a 2D U(I) QSL with spinon 
Fermi surfaces would show up. That is, the specific heat 
Cy ^ 7^2/3^ ^ Pauli-like spin susceptibility in the 

low temperature limifP. The crossover in the behaviour 
of the spin susceptibility from the local moment Curie- 
Weiss regime to the Pauli-like behaviour is expected to 
happen at the temperature s et by the bandwidth of ac¬ 
tive spinon bands (see Sec. IIIB), or equivalently, by 
the low-temperature Curie-Weiss temperature below T*. 
This crossover temperature should be very low because 
of the suppressed Curie-Weiss temperature at low tem¬ 
peratures. As a result, the Pauli-like spin susceptibility 
may be smeared out by various extrinsic factors like local 
magnetic impurities at very low temperatures. Likewise, 
even though the CyjT experiences a upturn below I OK 
in the absence of external magnetic fields, it is likely that 
the nuclear Schottky anomalies may complicate the spe¬ 
cific heat data. 


On the other hand, the apparently gapless spectrum 
of the spin excitations in the inelastic neutron scatter¬ 
ing measuremenfP3 is certainly consistent with the gap¬ 
less spinon Fermi surface of our U(I) QSL. Moreover, 
the measurements of relaxation rate from both NMR 
(I/(TiT)) and /rSR (AT“^) also indicate gapless spin- 
spin correlational^. In our U(I) QSL, the reduction of 
the spinon bandwidth due to the PCO increases the den¬ 
sity of the low-energy magnetic excitations. This would 
lead to a low-temperature upturn of the spin-lattice re¬ 
laxation, which is in fact observed in NMR and /iSR 
experimental^. 

A direct measurement of the PCO at low tempera¬ 
tures is crucial for our theory. To this end, a high resolu¬ 
tion X-ray scattering measurement and NQR (nuclear 
quadrupole resonance) can be helpful. Moreover, the 
presence of local quantum entanglement within the res¬ 
onant hexagon may be probed optically by measuring 
the local exciton excitations. Furthermore, if the sys¬ 
tem is in a U(I) QSL with a spinon Fermi surface, the 
low-temperature thermal conductivity can be an indirect 
probe of the low-energy spinon excitation, and a direct 
measurement of the correlation of the emergent U(I) 
gauge field might be possible because the strong spin- 
orbit coupling of the Mo atoms can enhance the coupling 
between the spin moment and the spin textur^^. 

A previous work on LiZn2Mo308 has proposed a the¬ 
ory based on varying spin exchange interaction from the 
emergent lattice that is caused by the lattice distortiorP^. 
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[Mo-Mo] u 

[Mo-Mo] d 

A e 

■/M 03 Ref 

LiZn2Mo308 

2 . 6 A 

3 . 2 A 

1.23 

7 

m 

Li2lnMo308 

2 . 54 A 

3.25 A 

1.28 

7 

139] 

ScZnMo308 

2.58A 

3.28A 

1.27 

7 

m 


TABLE L Mo-Mo bond lengths, anisotropic parameters (A) 
and number of electrons per M 03 O 13 cluster for three differ¬ 
ent cluster magnets. The electron number is counted from 
stoichiometry. 

In contrast, our work here is based on the electron 
degrees of freedom and their interactions. In Ap¬ 
pendix. we clarify the underlying quantum chemistry 
of LiZn 2 Mo 308 and justify the extended Hubbard model 

ofEq. Q. 

B. Other Mo based cluster magnets 



FIG. 7. (Color online.) (a) The M 03 O 13 cluster (adapted 
from Ref. 1151) . (b) The schematic energy level diagram of the 
molecular orbitals for a single M 03 O 13 cluster. The molecular 
orbitals are classified according to the irreducible representa¬ 
tions of the C 3 ^; point group of the clustei^. The unfilled 
molecular orbitals at high energies are not shown. 


The compounds that incorporate the M 03 O 13 cluster 
unit represent a new class of magnetic materials called 
“cluster magnets”. Several families of materials, such as 
M 2 M 03 O 8 (M = Mg, Mn, Fe, Co, Ni, Zn, Cd), LiRMo 308 
(R = Sc, Y, In, Sm, Gd, Tb, Dy, Ho, Er, Yb) and other 
related variants^^^®, fall into this class. The magnetic 
properties of most materials have not been carefully stud¬ 
ied so far. In Tab. |lj we list three cluster magnets with 
odd number of electrons in the M 03 O 13 cluster unit. We 
introduce a phenomenological parameter A to character¬ 
ize the anisotropy of the Mo Kagome lattice, which is de¬ 
fined as the ratio between inter-cluster (or down-triangle) 
and intra-cluster (or up-triangle) Mo-Mo bond lengths. 


[Mo-Mo] ^ 
[Mo-Mo] 


(28) 


Large anisotropy tends to reduce the interaction V 2 and 
increase the hopping ti so that the systems are more 
likely to be in the regime where the electron is only lo¬ 
calized on the up-triangle while the electron number on 
the down triangle is strongly fluctuating (see Fig.[^). In 
such a regime, there is no PCO, each localized electron on 
the up-triangle forms a local spin-1/2 moment, and these 
local spin-1/2 moments form a triangular lattice. If the 
system is in the weak Mott regime like the organic^^m^, 
the spin ground state is expected to be the U(l) QSL 
with a spinon Fermi surface^. 

As shown in Tab. |lj Li2lnMo308 has a larger 
anisotropic parameter than LiZn2Mo308. Unlike 
LiZn2Mo308 that has two Curie-Weiss regimes, the spin 
susceptibility of Li 2 lnMo 308 is instead characterized by 
one Curie-Weiss temperature ©cw = —207K down to 
25Kpni Moreover, the Curie constant is consistent with 
one local spin-1/2 moment in each up-triangle. Below 
25K, the spin susceptibility of Li2lnMo308 saturates to a 
constant, which is consistent with the Pauli-like spin sus¬ 
ceptibility for a spinon Fermi surface U(l) QSL. Besides 
the structural and spin susceptibility data, very little is 


known about Li2lnMo308. Thus, more experiments are 
needed to confirm the absence of magnetic ordering in 
Li 2 lnMo 308 and also to explore the magnetic properties 
of ScZnMo 308 and other cluster magnets. 
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Appendix A: Molecular orbitals and the Hubbard 
model 

As suggested by Refs. |T5] and [43l the Mo electrons 
in an isolated M 03 O 13 cluster form molecular orbitals 
because of the strong Mo-Mo bonding. Among the 7 
valence electrons in the cluster, 6 of them fill the lowest 
three molecular orbitals {A 2 , E^^^} in pairs, and the 

seventh electron remains unpaired in a totally symmetric 
Ai molecular orbital with equal contributions from all 
three Mo atoms (see Eig. [^. 

We first consider the molecular orbital states in the 
group {Ai, E^^\ This group can be described 

by a linear combination of an atomic state l^/i) at each 
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Mo site (which is in turn a linear combination of five Ad 
atomic orbitals). 

|Ai) = ^ [|' 0 i)a + |' 0 i)b + |' 0 i)c], (Al) 

lE^) = ^ [|^i)a + |^i)b + |V’i)c], (A2) 

|Ef) = ^ [|^i)a + |^i)b + |V’i)c], (A3) 

where /i(= A,B,C) labels the three Mo sites in the cluster 
and the atomic state is the contribution from the 

Mo atom at ja. The atomic states at different Mo 

sites are related by the 3 -fold rotation about the center of 
the cluster. Likewise, the fully-filled {A2, E^^^} and 

other unfilled molecular orbitals at higher energies are 
constructed from the atomic state |'02) and other atomic 
states \2pj) {j = 3 , 4 , 5 ), respectively. Here, the atomic 
states {\'ipj)f^} {j = 1 , 2 , 3 , 4 , 5 ) represent a distinct or¬ 
thonormal basis from the five Ad atomic orbitals that are 
the eigenstates of the local Hamiltonian of the MoOe oc¬ 
tahedron. 

We group the molecular orbitals based on the atomic 
state from which they are constructed. In this classifi¬ 
cation, for example, {Ai, E^^\ E^^^} fall into one group 

while {A2, E2^\ E2^^} fall into another group as they are 
constructed from two different atomic states. 

In LiZn2Mo308, the different molecular orbitals of the 
neighboring clusters M03O13 overlap and form molecu¬ 
lar bands. To understand how the molecular orbitals 
overlap with each other, we consider the wavefunction 
overlap of different atomic states Since the down- 

triangle has the same point group symmetry as the up- 
triangle in LiZn2Mo308, the wavefunction overlap of the 
atomic states in the down-triangles should approximately 
resemble the one in the up-triangles. More precisely, the 
wavefunction of the atomic state (e.g. IV^i)) has similar 
lobe orientations both inward into and outward from the 
M03O13 cluster, with different spatial extensions due to 
the asymmetry between up-triangles and down-triangles. 
Consequently, the orbital overlap between the molecular 
orbitals from the same group is much larger than the one 
between the molecular orbitals from the different groups. 
Therefore, each molecular band cannot be formed by one 
single molecular orbital but is always a strong mixture of 
the three molecular orbitals in the same group. 

We now single out the three molecular bands that are 
primarily formed by the group of {Ai, E^^^} molecu¬ 
lar orbitals. There are four energy scales associated with 
these three molecular orbitals and bands: 

( 1 ) the energy separation AE between the {Ai, 

(2') 

E\ } group and other groups of orbitals (both filled and 
unfilled), 

(2) the total bandwidth W of the three molecular bands 
formed by the {Ai, E^^\ E^^^} molecular orbitals, 

( 3 ) the intra-group interaction between two electrons on 
any one or two orbitals of the {Ai, E^^\ E^^} group. 


(4) the inter-group interaction between the electron on an 
orbital of the {Ai, E^^\ E^^^} group and the other elec¬ 
tron on an orbital of a different group. It is expected, 
from the previous wavefunction overlap argument, that 
the inter-group interaction is much weaker than the intra¬ 
group interaction and thus can be neglected at the first 
level of approximation. 

In this paper, we assume that the energy separation 
AE is larger than the total bandwidth W and the intra¬ 
group interaction. In this regime, the large AE separates 
these three molecular bands from other molecular bands 
(both filled and unfilled) so that the fully filled {A 2 , 
E 2 ^\ E 2 ^^} orbitals remain fully-filled and the unfilled 
molecular orbitals remain unfilled even after they form 
bands. Moreover, the large AE also prevents a band¬ 
filling reconstruction due to the interaction (in principle, 
the system can gain interaction energy by distributing 
the electrons evenly among different groups of orbitals). 
Therefore, we can ignore both the fully-filled and un¬ 
filled molecular bands and just focus on the three par¬ 
tially filled bands. It also means one will have to con¬ 
sider three-band model with all of {Ai, E^^} or¬ 

bitals on the triangular lattice formed by the M 03 O 13 
clusters. In this case, alternatively one could simply con¬ 
sider atomic states as the starting point. Then the rele¬ 
vant model would be a single-band Hubbard model based 
on the atomic state l^pi) at each Mo site of the anisotropic 
Kagome lattice. We take the latter approach in this pa¬ 
per. Einally, since only one atomic state \ 2 pi) is involved 
at each Mo site, the orbital angular momentum of the 
electrons are trivially quenched so that we can neglect 
the atomic spin-orbit coupling at the leading ordei®. 

The corresponding single-band Hubbard model is given 
by Eq. Q, where we include the on-site and nearest- 
neighbor electron interactions. Now it is clear that the 
physical meaning of the electron operator cj^ (^ia) 

Eq. 0 is to create (annihilate) an electron on the state 
with spin a at the Kagome lattice site i. 


Appendix B: Local moments in the strong PCO limit 


In Sec. we show the PCO reconstructs the spinon 
band structure and provide a possible explanation of 
the doule Curie regimes and 1/3 spin susceptibility in 
LiZn 2 Mo 308 . As we point out that the doule Curie 
regimes and 1/3 spin susceptibility are finite tempera¬ 
ture properties and independent from whether the spin 
ground state has a spinon band or not, in this section, 
we onsider an alternative and complementary strong cou¬ 
pling regime in which the PCO is strong such that the 
3 resonating electrons are almost fully localized in the 
resonating hexagon and form the local moments which 
then interact with each other. As far as the local mo¬ 
ment physics is concerned, the regime considered here is 
equivalent to the intermediate PCO regime in Sec. L^3t^. 

To elucidate the nature of the local moments in each 
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resonating hexagon, it is sufficient to isolate a single res¬ 
onating hexagon and understand the local quantum en¬ 
tanglement among the 3 resonating electrons. In this 
subsection, we consider two local interactions on the 
hexagon. The first interaction is already given in Eq. ^ 
which is the electron ring hopping model. The second 
interaction is the AFM exchange interaction between the 
electron spins. Since the electrons are always separated 
from each other by one lattice site, the AFM exchange is 
between the next nearest neighbors in the hexagon, 

= (Bl) 

im 

where are the lattice sites that refer to the 6 vertices 
of the resonating hexagon (see Fig. [^, rii is the electron 
occupation number at the site i and is the spin- 1/2 
operator of the electron spin at the site i. Because the 
electron position is not fixed in the hexagon, the AFM 
interaction is active only when both relevant sites are 
occupied by the electrons and we need to include rii into 
the exchange interaction. The full Hamiltonian for an 
individual resonating hexagon plaquette is composed of 
the above two interactions, 

Hp = -Y. [^l(4aC6a4/3C4/347C27 + 

a/37 

+IK2(cl„C2„4^C4^4^Ce^ + h.c.)] + (B2) 

Based on the perturbative values Ki = 6 tf/F 2 ^,IK 2 = 
Qt^lVi and the fact that ti > t 2 and Vi > V 2 in 
LiZn 2 Mo 308 , we think the relevant regime should be 
Ki ^ IK 2 . 

Because a strong PCO causes a strong modulation in 
the bond energy, the values of Ki and IK 2 for the resonat¬ 
ing hexagon would be modified from the effective Hamil¬ 
tonian that is obtained from the perturbative analysis. 
Likewise, the spin exchange in the resonating hexagon 
is enhanced from its perturbative value. So we expect 
this treatment is a good approximation to understand 
the local spin physics. 

The Hilbert space of the Hamiltonian is spanned 
by the electron states that are labelled by the posi¬ 
tions and the spins of the 3 resonating electrons. Be¬ 
cause the electrons are separated from each other by 
one lattice site, the Hilbert space for the positions is 
quite limited. There are a total of 16 states labelled 
by {|a/37)A = |2a,4^,67), |a/37)B = jla, 3/3, 57 )} with 
a, /3 ,7 =t, Since the local Hamiltonian commutes 
with the total electron spin Stot and 

to label the states. For both A and B elec¬ 
tron configurations, from the spin composition rule for 
3 spins (8) ^ (8) ^ | 0 | 0 |), we have 2 pairs of 

Stot = 1/2 states. 

The states with 5'tot = 3/2 are not favored by the 
AFM exchange. We find that when J > | (Ki 0 IK 2 — 

yKf - IK 1 IK 2 + K|) = IK 2 - ^ + the local 

ground states are 4 symmetric states with S'tot = 1/2. 



FIG. 8. (Color online.) Three singlet positions that are re¬ 
lated by the 3-fold rotation. 


Here, symmetric state refers to the symmetric linear su¬ 
perposition of spin states in A and B configurations. This 
is because the 3-electron Ki and IK 2 hopping terms hy¬ 
bridize the A and B configurations and favor symmetric 
superposition rather than antisymmetric superposition. 
The local 4-fold degeneracy can be effectively character¬ 
ized by 2 quantum numbers (s^,t^) with = ±1 and 
= ± 1 . refers to the total spin = 5'4»t ~ =*=0 
We also introduce a pseudospin- 1/2 operator r whose 
physical meaning is explained below. The wavefunctions 
of the four Ir^s^) states are given by (to the order of 
0 (K 2 /Ki)), 

in) = I [im)A - itn)A + int)B - im)B] (bs) 

lit) = ^ [2int)A - itn)A - im)A 

+2|tn)B - im)B - int)B] (b4) 

and Iti), \ii) are simply obtained by a time-reversal op¬ 
eration. 

What is the physical origin of this local 4-fold degener¬ 
acy? Clearly, the 2 -fold degeneracy of = ± 1/2 arises 
from the time-reversal symmetry and the Kramers’ the¬ 
orem. The remaining 2-fold degeneracy comes from the 
point group symmetry of the resonating hexagon. This 
is easy to see if we freeze the positions of the 3 electrons. 
To be concrete, let us fix the electrons to the sites 1,3,5 
in Fig. To optimize the exchange interaction, 2 elec¬ 
trons must form a spin singlet, which leaves the remain¬ 
ing electron as a dangling spin- 1/2 moment. As shown 
in Fig. this singlet can be formed between any pair 
of the electrons and the different locations of the spin 
singlet are related by the 3-fold rotation. Even though 
there seems to be 3 possible singlet positions, only 2 of 
them are linearly independent, which gives to the 2 -fold 
degeneracy which survives even when the ring electron 
hopping is turned on. As a result, the pseudospin r is 
even under time-reversal and acts on the space of the sin¬ 
glet position or equivalently the dangling spin position. 
In fact, the two states in Eqs. ( |B3[ ) and ( |B4[ ) comprise 
the E irreducible representation of the Cs^ point group. 

Since the spin s (pseudospin r) is odd (even) under 
time reversal symmetry, when the external magnetic field 
is applied to the system, only the spin s (s = 1 / 2 ) couples 
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to the magnetic field. Therefore, the three electrons in where we have defined the x variable as 

the resonating hexagon behave like one spin s = 1/2 in 

the magnetic field. This is how the Curie constant with 

the PCO becomes 1/3 of the Curie constant without the 

PCO. 


Appendix C: Levin-Wen’s variational dimer 
wavefunction approach 


= (C4) 

a 


Here we explain the string mean-field theory in 
sec. inra in details. To solve the combined Hamiltonian 
of Hs and ^ring self-consistently, we obtain the effective 
spinon hoppings by evaluating the boson or rotor hopping 
amplitudes with respect to the variational dimer ground 
state 

(l+(R)l;(r)) = (l+(R))(l;(r)) (Ci) 

where /r, label the sublattices. We also evaluate the 
parameter Mijkimn against the spinon hopping Hamilto¬ 
nian. Using the Wick theorem, we have 


and we have also used the three-fold rotational symme¬ 
try as well as the reflection symmetry of the hexagon in 
Eq.(|C3|) so that 


Mijkimn — 

a/37 

a 

+ J2^fLfla)ifmJja)ifljna) 

a 

cxf3 


al3 


-T.(fLfna){fL0fj0){fyi0) (C 2 ) 

al3 


Xij + ^ ^ “ -^XijXikXiU 


(C3) 


Xijf Xkl Xmn 

Xik = Xkm = Xmi = Xjl = Xln = Xnj (C6) 
Xik — Xjm — Xkn‘ 
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